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ABSTRACT 

^ Detection of the tensor perturbations predicted in inflationary models is im- 

f- — portant for testing inflation as well as for reconstructing the inflationary 

O potential. We show that because of cosmic variance the tensor contribution 

to the square of the CBR quadrupole anisotropy must be greater than about 

20% of the scalar contribution to ensure a statistically significant detection 

Q^ of tensor perturbations. This sensitivity could be achieved by full-sky mea- 

Q surements on angular scales of 3° and 0.5°. 



Inflation has had almost as much impact on cosmology as the big-bang 
model itself. It provides an explanation for the flatness and smoothness of the 
observed Universe, as well as the inhomogeneity needed to seed the structure 
seen today. However, inflation has yet to be tested in a signiflcant way. The 
keys to doing so are its three robust predictions: spatially-flat Universe (total 
energy density, including matter, radiation and vacuum energy equal to the 
critical density) and nearly scale-invariant spectra of scalar (density) [1] and 
tensor (gravity-wave) [2] metric perturbations. There is some evidence that 
the density of matter is close to unity [3] and several large-scale experiments 
are underway to search for the nonbaryonic dark matter that must exist if 
f^matter = 1- Scalar perturbations seed the formation of structure, and hence, 
observational data concerning the distribution of matter in the Universe can 
provide information about them. Today, tensor perturbations correspond to 
a stochastic background of gravitational waves that could possibly be de- 
tected by space-based interferometers. Both tensor and scalar perturbations 
give rise to anisotropy in the temperature of the Cosmic Background Ra- 
diation (CBR) on angular scales from arcminutes to 180° (see Fig. 1), and 
CBR-anisotropy is a promising means of detecting them [4]. 

Long before inflation the attractiveness of scale-invariant perturbations 
and a flat Universe had been emphasized [5]; thus tensor perturbations play 
a crucial role in discriminating between inflation and other "attractive" the- 
ories. Further, detection of the tensor perturbations is vital to the recon- 
struction of the inflationary potential [6]. Denoting the spectral indices of 
the scalar- and tensor-metric perturbations by n and rix respectively (scale 
invariance corresponds to n — 1 = nx = 0), and their contributions to the 
variance of the quadrupole CBR anisotropy by S and T, the value of the 
inflationary potential and its flrst two derivatives are given by [7] 



V{(f)5o)/mpi'^ = 


= 1.65(1 -1.4nT)r, 


U'(^50)/?7^P1^ = 


= ±S.3f-T. 


V"{(l)5o)/mpi^ = 


= 21[(n-l) + 0.4r] T 
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where r = -g. Subscript '50' indicates the value of the inflaton fleld (f) when 
present horizon-sized metric fluctuations crossed outside the horizon during 
inflation; this typically occurs when the scale factor was e^° smaller than its 
value at the end of inflation. In addition, there is a consistency relation: 
riT = —kr, which is an important test of inflation [8]. 
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CBR-anisotropy measurements may offer the best possibifity for reveafing 
tlie presence of tensor perturbations [9]. Because tlie predictions for tlie 
metric perturbations are statistical in nature and tlie sky is but a finite 
sample of the Universe, sampling variance, or cosmic variance [10], provides 
a fundamental limit to the separation of the tensor and scalar contribution 
to CBR anisotropy. In this Letter we show that sampling variance precludes 
the detection of tensor perturbations if r < 0.2, and further, we show how 
a simple experiment involving measurements of CBR anisotropy on angular 
scales of around 3° and 0.5° can achieve this limiting sensitivity. 

At present, imprecise knowledge of cosmological parameters, Hubble con- 
stant. Ho = lOO/ikms"^ Mpc~ , baryon fraction, 0^, scalar spectral index n, 
and cosmological constant, as well as the ionization history of the Universe 
limit the sensitivity further, as has been emphasized in Refs. [11]. Further, 
anisotropy measurements are still dominated by experimental error rather 
than cosmic variance. We believe that prospects for determining the cos- 
mological parameters by other means is good and experiments are rapidly 
becoming more accurate, so we take the optimistic view that the limiting 
factor may soon be sampling variance. For definiteness, we take h = 0.5, 
fts = 0.05, A = 0, n = 1, and assume standard ionization history. 

We begin by briefly reviewing the statistics of CBR anisotropy in gen- 
eral [12], and that produced by inflation in particular. CBR-temperature 
fluctuations are expanded in spherical harmonics, 

6T{e,<P) = J2ai^YUe,<P). (2) 

Isotropy in the mean guarantees that (a/m«*'m') — CiSniSmm'^ where brackets 
indicate average over an ensemble of observers. It is the variance of the 
multipoles that encodes information about the metric perturbations, and 
Ci = (|a/mP) is called the angular-power spectrum. (The expectation for 
the square of the quadrupole anisotropy Q^ = 56*2 /47r.) Provided that the 
underlying perturbations are Gaussian (as is almost certainly the case for 
inflation), all predictions can be derived from the angular-power spectrum. 
With access to only one sky, for each / we can only measure 2/-|-l independent 
multipoles. Thus, the estimator for C;, Cf ^ = J2^ |a;mP/(2/ -|- 1), will differ 
from Ci due to finite sampling: 

((Cf ^ - Cif) = ^^^^^^^^ = J^YJl^- (3) 



This ultimate uncertainty in our knowledge of the angular-power spectrum 
is why there is a limiting sensitivity to r. 

Most experiments do not directly measure the angular-power spectrum, 
but instead, the square temperature fluctuation on a given angular scale, 
whose expectation is 

{ST^)=J:^CW. (4) 

The window function Wi depends on the beam size and chopping strategy. 
Very roughly, an experiment that measures the temperature difference be- 
tween directions separated by angle 6 with beam size cr ~ ^ has a window 
function that is centered around / ~ vr/^, with width of order /. If the 
experiment samples the full sky, then the variance of ((^r^)sky is 

{m%,,-{6T^)r) = 2j:^^cfw^. (5) 

The subscript "sky" distinguishes sky average from ensemble average {ST'^). 

Figure I shows the angular-power spectra arising from scale-invariant 
scalar and tensor perturbations [13]. The expectation for an experiment is 
determined by the sum of the scalar and tensor contributions, Ci = Cf -\-Ci ■ 
(Likewise, Q"^ = S -\- T, where S = 5Cf /47r and T = 5Cj /Air.) Differences 
in the metric between two points on the last-scattering surface (Sachs- Wolfe 
effect) are the dominant contribution to both angular-power spectra at large 
angles (/ < 60); tensor perturbations decrease in amplitude once their wave- 
lengths become smaller than the Hubble radius which explains why the tensor 
spectrum drops at small angles (/ > 60) [14]. The scalar angular-power spec- 
trum increases slowly with / due to the pre-recombination oscillations of the 
baryon-photon fluid, giving rise to the "Doppler peak" at / ~ 200 [15]. 

A simple strategy for separating Cf and Cf is to measure temperature 
fluctuations at two angular scales. The window function for one experiment, 
Wi^ , is centered at small angles (/ ~ 200) where Cf is insignificant. The 
window function for the other experiment, Wj^, is centered at intermediate 
angles (/ ~ 50). In essence, experiment B measures the amplitude of the 
scalar spectrum, and experiment A measures "excess power" at small / due 
to tensor perturbations. We define a measure of this excess power. 
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Figure 1: Angular-power spectra arising from scalar and tensor metric per- 
turbations for n — 1 = nx = 0, /i = 0.5, and fts = 0.05. The dashed curve is 
the optimal window function, Wj^, described in the text. 



where ((5r|)r=o = J2i 'T~^T^i i^ the anisotropy expected in experiment A 
for r = and (<^r|)sky is that observed. Note that Z is proportional to r. 

Because of cosmic variance a measurement of Z ^ does not necessarily 
imply r ^ 0. Consider the cosmic variance in Z: 



AZ' _ {{Z-{Z)f) 2Y.i[{2l + l){Cf{r) + Crf{W{'f] 



{zy 



{zy 



{T.im+l)Cf{r)Wi']f 
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(We have not included the cosmic variance from experiment B because in the 
case of a full-sky measurement it is negligible.) The ratio Z/AZ, which is 
proportional to r, is a measure of signal-to-noise. Choosing a standard form 
for the window function, W{^ = [1 ~ P;(cos ^)] exp( — /(/ -|- l)cr^), which cor- 
responds to Gaussian beamwidth a and chop angle 6, we maximized ZjAZ 
(i.e., sensitivity to r) by varying Q and a . The optimal window peaks around 
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Figure 2: Probability distribution P{Z\r) for (n, r) = (1, 0), (1, 0.1), (1, 0.2) 
and (1, 0.3) (solid); for (n, r) = (0.94, 0) (broken); and for (n, r, fts) = (1, 
0, 0.045) (dot-dash). 



/ = 50; see Fig. 1. Its shape can be understood in terms of signal-to-noise 
weighting: For small / there is large noise (cosmic variance); for large / there 
is small signal {Cf jC^' is small). 

We have numerically calculated the probability distribution P[Z\r) using 
a Monte-Carlo method; see Fig. 2. For r ~ 0.3 it is clearly possible to rule 
out the hypothesis that r = 0, since its probability distribution has little 
overlap with that for r = 0.3. 

To address the overlap of distributions quantitatively, statisticians define 
"size" and "power". The power is the probability, of measuring Z < Z\ 
given r = 0, and the size is the probability of measuring Z < Z\ given 
r = r\. Figure 3 illustrates the likelihood for ruling out r = with 95 per 
cent confidence, given that the actual value r 7^ 0. In other words, we have 
fixed the power to be 0.95, and calculated the size as a function of r. As r 
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Figure 3: One minus the fraction of sky realizations (= size) that lead to 
a 95 per cent confidence (power = 0.95) or 99.7 per cent confidence (power 
= 0.997) exclusion of the r = hypothesis as a function of r. The broken 
curves are for the Z-test, and the solid curves are for the likelihood-ratio 



test; n 



1 (open symbols) and n = 0.85 (filled symbols). 



increases, P[Z\r) overlaps less with the distribution for r = 0, and thus the 
size decreases. For r ~ 0.2 there is a 95 per cent chance that one can exclude 
the r = hypothesis with at least 95 per cent confidence. For the infiationist 
who feels lucky, for r = 0.1 there is a 50 per cent chance of being able to 
eliminate the r = hypothesis with at least 95 per cent confidence. 

While we have chosen Wj^ to maximize Z/ AZ , the breadth of this maxi- 
mum (in window-function space) is large so that even substantial changes in 
W{^ do not greatly affect Z/ AZ . One might ask if a more cleverly defined 
observable could achieve a greater sensitivity. To the contrary, we now show 
that the sensitivity cannot be significantly improved. 

To address this issue we apply a likelihood-ratio test [16] that uses all 
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possible information, i.e., all the multipoles. Likelihood-ratio tests are de- 
signed to discriminate between two hypotheses, Hq and Hi. Their virtue 
is that they are "most powerful;" that is, for fixed size such a test yields 
maximum power. Let hypothesis Hq be the assertion that r = 0, and Hi be 
the assertion that r = ri. 

We define the likelihood ratio, 

, _ P{{aim]\r = ri) __,.^( Cf V"^'^^' 



A* = E 



P({a,„}|r = 0) t\\Ct + Cr{ri), 

E|„ 12 / 1 1 

m \'^lm\ /J- J- 



1=2 




c?+cr{ri) 



(9) 



We then proceed as before: generation of the distributions P[X\r = 0) and 
P[X\r = ri), followed by calculation of size for a given power. The results are 
shown in Fig. 3 and are nearly identical to those of the previous "Z-test." 

In principle, all the multipoles relevant for the likelihood-ratio test could 
be obtained from an all-sky experiment on the 0.5° scale. This does not 
mean that the virtues of the "Z-test" are solely pedagogical. Its sensitivity 
is almost the same as the likelihood-ratio test, and full-sky coverage may 
not be possible, making the extraction of multipoles from any experiment 
difficult. In the case of partial-sky coverage, the "Z-test" may offer a practical 
alternative whose limiting sensitivity scales roughly as [1//a + 1/(16/b)]^'^ 
where /a and fs are the fractional sky coverages of experiments A and B. 

Depending on the actual values of /i, 0^, A, and n the ultimate sensitivity 
to r varies; it is most sensitive to the shape of the scalar spectrum on the 
Sachs- Wolfe plateau — and hence n. Decreasing n from 1 to 0.85 improves 
the limiting sensitivity to r ~ 0.1 (see Fig. 3); increasing n has the opposite 
effect. The limiting sensitivity is insensitive to fts and h: for fts = 0.02 it 
decreases to r ~ 0.19, and for h = 0.4 it increases to r ~ 0.21. 

What about the effect of imprecise knowledge of cosmological parameters? 
Figure 2 shows that it is not possible to distinguish n = 0.94, r = from 
n = 1.0, r = 0.2. Suppose n ~ 1, but one could not preclude a value as 
small as 0.94, the limiting sensitivity increases to r ~ 0.4. Inferring extra 
power on the Sachs- Wolfe plateau (i.e., r ^ 0) requires accurate knowledge 
of the height of the Doppler peak, which depends upon fts and h. Thus, 
uncertainty in fts and h can mimic r ^ 0. For example, fts = 0.045, 
r = mimics fts = 0.05, r = 0.04 (see Fig. 2), and a 10% uncertainty in 
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fts increases the limiting sensitivity to about 0.24. The effect for h is even 
greater: h = 0.55, r = cannot be distinguished from h = 0.5, r = 0.f4. 

A logical extension of the two-scale measurement is the extraction of 
both r and rix via a three-scale measurement. An independent measure of 
riT permits the testing of the consistency relation: r = —Tut- However, even 
for r as large as 1, so that rix = —1/7, ruling out nx = is nearly impossible. 
For r < 1, falsification of the consistency relation requires |nx| ^ r/7 [17]. 

Finally, while we have focussed solely upon anisotropy, polarization is 
another way of separating tensor and scalar perturbations. For / < 10, 
linear polarization of CBR anisotropy due to tensor perturbations is about 
\/60r times greater than that due to scalar perturbations [18]. By means 
of a likelihood-ratio test for polarization and assuming an experiment only 
limited by cosmic variance, we find a sensitivity of r ~ 0.03. However, the 
polarization-to-anisotropy ratio is less than 1 per cent for r = 1 — and closer 
to 0.01% for r ~ 0.03 — and polarization is difficult to measure on large- 
angular scales. For the foreseeable future it seems unlikely that polarization 
experiments can achieve higher sensitivity. 

In summary, CBR anisotropy appears to be the most promising means 
of discovering infiation-produced tensor perturbations. We have shown that 
even with perfect knowledge of the temperature of our CBR sky and the 
cosmological parameters, the ratio of tensor to scalar perturbations r must be 
greater than about 0.2 to guarantee a statistically significant detection. The 
importance of discovering tensor perturbations provides further motivation 
for the cosmic-variance limited, full-sky CBR anisotropy experiments being 
planned for satellite platforms. 
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